d-wave resonating valence bond states of fermionic atoms in optical lattices 



O 
O 



Simon Trebst 1 ' 2,3 , Ulrich Schollwock 4 , Matthias Troyer 1 , Peter Zoller 5 
^ Theoretische Physik, Eidgenossische Technische Hochschule Zurich, CH-8093 Zurich, Switzerland 
( 2 ' Computational Laboratory, Eidgenossische Technische Hochschule Zurich, CH-8092 Zurich, Switzerland 
Microsoft Research and Kavli Institute for Theoretical Physics, 
University of California, Santa Barbara, CA 93106, USA 
^ Institut fur Theoretische Physik C, RWTH Aachen, D-52056 Aachen, Germany and 
( 5 ' Institute for Theoretical Physics, University of Innsbruck, and 
Institute for Quantum Optics and Quantum Information of the Austrian Academy of Science, 6020 Innsbruck, Austria 

(Dated: February 2, 2008) 

We study controlled generation and measurement of superfluid d-wave resonating valence bond 
(RVB) states of fermionic atoms in 2D optical lattices. Starting from loading spatial and spin 
patterns of atoms in optical superlattices as pure quantum states from a Fermi gas, we adiabatically 
transform this state to an RVB state by change of the lattice parameters. Results of exact time- 
dependent numerical studies for ladders systems are presented, suggesting generation of RVB states 
on timescale smaller than typical experimental decoherence times. 



CO 

c3 



i 

C 

O 
o 



(N 
> 

o 

oo 

o 
in 
o 

c3 



o 
o 



X 



Resonating valence bond (RVB) states, in which elec- 
trons are paired into short-range singlets by a purely 
electronic mechanism, were originally introduced [1| as 
wave functions for insulating spin liquid ground states 
of Mott insulators. Shortly after the observation of 
high- temperature superconductivity in the cuprates 0, 
Anderson conjectured that they might be described by 
doped RVB states Q, in which the RVB pairs of the in- 
sulating state become mobile and superconducting upon 
doping. Since the predicted exotic d-wave pairing sym- 
metry of the RVB pairs was confirmed experimentally 
in the cuprate superconductors p| , the RVB scenario re- 
mains as one of the most promising contenders for the 
theory of high temperature superconductivity Q . While 
the ground state of the lightly doped two-dimensional 
Hubbard model is still unknown, the d-wave RVB sce- 
nario has been confirmed for t-J and Hubbard models of 
coupled plaquettes 0, and ladders (consisting of two 
coupled chains) || by numerical simulations. The key 
question is whether the RVB state of weakly coupled lad- 
ders or plaquettes [8j is connected to a possible RVB 
ground state of the square lattice Hubbard model: does 
the RVB state survive when the inter-plaquette (ladder) 
coupling is increased and becomes the same as the intra- 
plaquette (ladder) coupling, at which point we have a 
uniform square lattice? 

We propose to address this question using cold atoms 
loaded in optical lattices, which allow the realization of 
Hubbard models with controllable parameters, promis- 
ing an entirely new avenue in the s tudy of strongly cor- 
related systems [HElElElElEiElEi- In par- 
ticular, formation of RVB ground states with fermionic 
atoms in a 2D geometry is one of the ultimate challenges 
|l2| . However, d-wave pair binding energies are two or- 
ders of magnitude smaller than the hopping amplitudes 
for atoms in an optical lattice which are typically in 
the kHz range: this poses a significant experimental chal- 
lenge in terms of temperature requirements etc., which 



will be difficult to meet. Instead we propose here the 
formation of RVB states by loading spatial and spin pat- 
terns of atoms in optical superlattices as pure quantum 
states from a reservoir of a quantum degenerate Fermi 
gas |l9f . and adiabatically transforming this state to an 
RVB state by change of the lattice parameters. An initial 
lattice configuration is designed in such a way that the 
desired initial state corresponds to the ground state of 
the deformed lattice in the form of a simple many-body 
product state, with a large excitation gap. This makes 
the preparation of these states robust against imperfec- 
tions. In choosing a protocol for the deformation of the 
lattice to achieve the RVB state, the challenge is to se- 
lect a parameter path, which minimizes the number of 
possible avoided crossings, and which sets the time scale 
for the adiabatic transformation, reminiscent of discus- 
sions in adiabatic quantum computing [2(J. While this 
technique is illustrated for the case of RVB states here, it 
might provide a generic procedure to generate non-trivial 
ground states of dilute fermionic gases in lattices. 

In this Letter we present a detailed investigation of this 
scenario. We will start with the loading of spin and spa- 
tial patterns on decoupled plaquettes formed by optical 
superlattices representing different atomic dopings, and 
discuss requirements for the formation of RVB plaquette 
states, and possible experimental signatures demonstrat- 
ing d-wave pairing. We present results of exact time- 
dependent numerical studies for ladders systems, sug- 
gesting that RVB states can be generated on a timescale 
smaller than the typical decoherence time of atoms in 
optical lattices. Since the numerical simulation of 2D 
strongly correlated systems belongs to the class of non- 
deterministic polynomial-time hard (NP-hard) computa- 
tional problems |2l|. we are not able to provide answers 
for the coupling of plaquettes or ladders to probe the 
ground state of the 2D Hubbard model, but pose this as 
an important problem left for experiments to be solved. 

Optical lattices - We start by summarizing the basic 
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FIG. 1: Optical potentials generated by two counterprop- 
agating laser beams with wave vector k — 2ir/\. a) V = 
V sm{kx) 2 ; b) V = Vbsin(fcx) 2 + Visin(fcx/2) 2 ; c) V = 
V sm(kx) 2 + Vi sm(kx/2) 2 + V 2 x; d) V = Vo sin(kx) 2 + V 2 x. 



tools available for construction of ID ladders and 2D 
square lattices with optical potentials [ill Il2l IT3L [iH Hif . 

Counterpropagating laser beams along the Xi (i = 1, 2, 3) 
directions realize a periodic optical potential Vq(x) — 
Voi sin(fca;i) 2 with wave vector fc = 2ir/\, optical wave 
length A, and lattice depths Vbi controlled by the laser in- 
tensities and atomic alternating current (AC) polarizabil- 
ities. These potentials are additive, provided the lasers 
generating the lattice have slightly different optical fre- 
quencies, which are easily generated as sidebands from 
the original laser beam, so that the interference terms 
average out. Furthermore, interference of two lasers with 
angles ±9 between their propagation directions give rise 
to potentials of the form Vu s'm(k'xi + <&) 2 with effec- 
tive wave vector k' — k cos 9 and displaced by a phase 
<f>. This allows us to construct in particular potentials 
Vi(x) = YJ i Vu sm(kxi/2) 2 . Adding these various poten- 
tials gives rise to superlattice structures. We can also add 
potentials in the form of linear ramps, V^ix) = VnXi 
by sitting in the wing of a focused laser beam, as in far-off 
resonance laser traps. Fig. ^ illustrates lattice configu- 
rations which can be generated in this way by varying 
the intensity parameters Vb,i,2i which will be employed 
below for the construction of d-wave pairing. 

The dynamics of cold atoms confined to these optical 
potentials is described by a Hubbard model 0, 0, 0, 

mm 

H = - ^ til (4a C 3 a + h - C )+ U Y n i 1 i n n+^2^ n ^ ' 

i i,<7 

where the c, lCT are fermionic annihilation operators with 
spin and n^o- = cj a c i a is the particle number on site i. 
The tij are spatially dependent hopping matrix elements 
connecting neighboring sites i and j, and the fii are site 
offsets, as determined by the superlattice structure (see 
Fig. QJ. The collisional repulsion U between the atoms 
can be controlled by Feshbach resonances [l^. E^ . In 
writing the Hubbard model we have assumed that the op- 
tical potentials are spin independent, which is the usual 
case of Alkali atoms in their ground state 0] ■ 

Plaquette RVB States - The motion of the atoms can 
be confined to a 2D lattice by a strong transverse opti- 
cal potential Vq z . Employing superlattices we can gen- 




FIG. 2: Schematic illustration of the adiabatic protocol that 
generates RVB states on two decoupled plaquettes with 4 
atoms (left) and 2 atoms (right). The protocol sequence is: 
1) t — > 1; 2) fj, — > 0; 3) fi± — » 0. Optical potentials are 
sketched for x and ^-directions and the large (intermediate, 
small) circles indicate two (one, zero) atoms on a site. 



erate double well potentials corresponding to decoupled 
plaquettes. Our first goal will be to study atomic dy- 
namics on these uncoupled plaquettes, and in particular 
the generation of RVB (ground) states. The strategy 
is to (i) deform the optical lattice on each plaquette so 
that the corresponding ground state has a simple product 
form, \ip) = rii<T C LlO)' an d ( n ) adiabatically transform 
the lattice into an unperturbed plaquette, so that the fi- 
nal ground state is the desired RVB state. To prepare 
such a pure initial state we load atoms from a reservoir 
of a quantum degenerate Fermi gas via a coherent or dis- 
sipative Raman process into the optical lattice [2^]. By 
choosing an appropriate pattern of site offsets in the lat- 
tice, we ensure only atoms are transferred which match 
the energy conservation condition. As shown in |2^|, this 
allows us to filter out from an initial finite entropy en- 
semble a pure spin and spatial pattern of atoms. 

We consider preparation of one plaquette with 4 atoms 
(half-filling) and one plaquette with 2 atoms. The initial 
optical lattice for these two cases is described in Fig.|2Ji), 
which amounts in the limit t± <C to a preparation of 
the states with 4 (2) atoms in a product state of the form 

14)^=^4^10) , 12)^=4^10) , (1) 

where the indices 1 and 3 denote the wells along the x- 
direction as illustrated in Fig. [21 

An adiabatic protocol that allows to transform this ini- 
tial product state to an RVB plaquette state is presented 
in Fig. [3 using the transformations of optical lattices il- 
lustrated in Fig. ^ Full plaquette symmetry is restored 
by coupling two wells first along the x-direction with 
their fully depleted counterparts along the y-direction. 
The initial optical lattice breaks reflection symmetry 
along the x-direction which upon coupling by increasing 
t (Vq x — > 0) enforces that the state |2)^ adiabatically 
connects only to the dimer state which is antisymmet- 
ric under spin exchange and symmetric under reflection. 
Eliminating the chemical potential shifts /1 (V2 X —> 0) and 
A*_l — > 0) one can adiabatically prepare the ground 
states of a 4-site plaquette state. For 4 atoms it is given 



3 




HQ 



EC 



128 256 384 

T [lit] 



B til 



FIG. 3: Signature for d-wave RVB pairing when decoupling 
two plaquettes with 6 atoms for varying onsite repulsion U ft. 
b) For U /t < 4.5 pairs are formed with a small binding energy 
shown in the inset. In the final state after sufficiently slow 
decoupling the hole pair is located on one plaquette, while for 
large repulsion U/t > 4.5 the unpaired holes separate into 3 
atoms on each plaquette (right panel), a) Projection of the 
final state onto the subspace with 3 atoms on each plaquette. 
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where we have only written the dominant terms omit- 
ting states with local double occupancy, and 
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While the state |2) has s-wave symmetry with respect 
to the vacuum, it has d x 2_ y 2 -symmetry with respect to 
the ground state |4) at half- filling, as the bosonic d-wave 
hole pair operator A ~ (si,2 + S3, 4 — si,4 — ^2,3) /2 has a 
non-vanishing contribution (2| A |4) ^ 0. The state |2) is 
therefore referred to as the d-wave RVB state (of the hole 
pair). The whole transformation is protected by a finite 
gap that sets the time scale for an adiabatic coupling 
with a fidelity larger than 0.99 to T w 50/i 

Coupling of plaquettes - Coupling and decoupling these 
two plaquettes can provide a clear experimental signa- 
ture for the pairing of atoms (or holes) in these RVB 
states, see Fig.|3] When coupling the plaquettes the holes 
move to the center of the 8 coupled sites as illustrated 
in Fig.|U>). Subsequently decoupling the two plaquettes, 
the various atomic occupation numbers on the two pla- 
quettes are in direct competition. If a bound state of 
two atoms (holes) is disfavored by large onsite repulsion 
(inset of Fig.|3^), the atoms will maximize their respec- 
tive kinetic energies as the plaquettes are decoupled and 
the holes separate. After adiabatic decoupling there are 



two plaquettes with 3 atoms each as shown in Fig. [3J If 
a bound state is favorable the pair hops to either one of 
the plaquettes, leaving the system in a final state with an 
even number of atoms on both plaquettes. Pair binding 
is found for U/t < 4.5 with a maximum bindingenergy 
of E b /t 0.04 for U/t » 2.5 (see inset of Fig. 01 8]. The 
time scale for the adiabatic decoupling is dominated by 
the small binding energy scale and found to be T sa 500/i 
or T « 1/2 seconds for 40 K atoms which is well within 
the decoherence time of these systems |l6j . Reversing 
the adiabatic protocol illustrated in Fig. [3 the plaquette 
states after decoupling are transformed into single-site 
states. For pair binding the system ends up with 3 dou- 
bly occupied sites, while for unbound pairs there will be 
two doubly occupied sites and two sites with one atom 
each. Experimentally the two scenarios can be distin- 
guished by associating the atoms into molecules 0, 
with the number of formed molecules having a ratio of 
3/2 respectively. 

The detailed adiabatic protocol that couples two pla- 
quettes is illustrated on top of Fig. ^ While the final 
ground state of two plaquettes with 6 atoms is symmet- 
ric with respect to the exchange of the two plaquettes, 
coupling two plaquettes prepared in the states |4) and |2) 
will adiabatically connect the initial state |4) |2) to both 
(anti)symmetric combinations (|4) |2) ± |2) |4))/v2- In 
order to prevent coupling to the antisymmetric state we 
need to explicitly break reflection symmetry along the in- 
direction. This can be achieved by ramping the chemical 
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FIG. 4: Singlet gaps for the path that couples and decou- 
ples two plaquettes illustrated by the optical potentials on 
top. Two plaquettes with 4 and 2 atoms are coupled for 
U/t = 2.5 with a chemical potential gradient fi being applied 
along the ^-direction (first three panels). The decoupling of 
plaquettes is shown for bound (fourth panel) and unbound 
pairs (right panels). Dashed (dotted) lines indicate states to 
which a transition from the ground state is forbidden due to 
symmetry constraints (particle conservation on the individual 
plaquettes). The x- and {/-parities are given by "+" (even) 
and "-" (odd). 
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FIG. 5: Ramping times needed to achieve a fidelity larger than 
0.999 for an adiabatically evolved wavefunction as the chemi- 
cal potential shift /i at predefined lattice sites is reduced. The 
adiabatic protocol is performed in small sequences /j, — » fi — S/j, 
with <5/i = 0.1. Doping S — 1/8. On the right: a) decoupled 
plaquettes; b) comb-like structure; c) ladder geometry. 

potential as illustrated in Fig.QJ;). If we subsequently in- 
crease the hopping between the two plaquettes (Vi x — > 0) 
and finally eliminate the shift in the chemical potential 
(V2x —> 0), we can adiabatically connect the initial state 
to the symmetric ground state of two coupled plaquettes. 
The coupling transformation is protected by a consider- 
able gap as shown in Fig. 0] The gaps for the subsequent 
decoupling process which keeps full exchange symmetry 
between the two plaquettes reveal the small pair bind- 
ing energy in the vicinity of small interplaquette hopping 
t as shown in Fig. 01 for U/t = 2.5 (pair binding) and 
U/t = 6.0 (no binding). The low-energy dynamics shown 
in Fig.0|can be probed by measuring the structure factor 
in light-scattering spectroscopy 25]. 

Doped d-wave RVB ladders - Finally, we discuss half- 
filled ladder systems by coupling multiple plaquettes and 
the preparation of a hole doped <i-wave RVB state, d- 
wave hole pairing in ladders is strongest for U/t w 2.5; 
pairs are localized and pinned with a period 1 / 8 for open 
boundary conditions and doping 8. To enable a fast adia- 
batic coupling process we prepare an initial optical lattice 
that already places hole pairs close to their final location 
in the doped ladder by mimicking the 1/8 periodicity as 
shown in Fig.Et)- Using a pattern loading technique 
we can load a gas of fermionic atoms into this optical lat- 
tice so that the shifted sites in this comb-like structure 
remain depleted while all other rungs in the ladder sys- 
tem are half-filled. When adiabatically coupling this ini- 
tial state to the final RVB state of the doped half-filled 
ladder by reducing the chemical potential shifts \x, the 
time scale is dominated by the time needed to establish 
phase coherence between previously unconnected ladder 
parts. 

In our numerical simulations, we considered a ladder 
with 32 rungs filled with 56 particles, i.e. doping (5 = 1/8. 
In the initial configuration 28 rungs are half-filled and 
the rungs 4, 12, 20, 28 depleted, see Fig. |3]d). Using 
adaptive time-dependent density matrix renormalization 



group (DMRG) algorithms (2(|, we have calculated the 
time scales for ramping down fj, at the specified sites. 
We find that with an overall ramping time T w 250/t, 
following the protocol shown in Fig. EJ, the fidelity is 
larger than 0.9 . To improve this bound to 0.99, we 
estimate T ~ 450/t. 

Phase coherence in the ladder RVB state can be ex- 
perimentally probed via molecule formation where the 
molecules are generated through a laser-induced Ra- 
man transition that couples to a molecular state m with 
d x 2_ y 2 -symmetry with a transition matrix element of the 
form m ' A jl5j . If the ladder exhibits quasi-long range or- 
der the generated molecules will be phase coherent form- 
ing a quasi-condensate. A truly long-range superfluid 
ground state of ci-wave RVB pairs can be stabilized in 
experiment by weak coupling of ladders. Increasing the 
inter-ladder coupling in the experimental setup will an- 
swer one of the foremost open questions of solid state 
physics: is the d-wave RVB state on ladders adiabatically 
connected to the ground state of the Hubbard model on 
the uniform square lattice, and hence the RVB theory 
of high temperature superconductivity confirmed - or is 
there a quantum phase transition to a new phase? 
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